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o ■ 

^ ■ ABSTRACT 

C/3 ', 

, ^ , ' Aims. We theoretically and phenomenologically investigate the question whether the 7-ray emission from the remnants of the type la super- 
novae SN 1006, Tycho's SN and Kepler's SN can be the result of electron acceleration alone. 
^-H • Methods. The observed synchrotron spectra of the three remnants are used to determine the average momentum distribution of nonther- 
^ ' mal electrons as a function of the assumed magnetic field strength. Then the inverse Compton emission spectrum in the Cosmic Microwave 
Background photon field is calculated and compared with the existing upper limits for the very high energy 7-ray flux from these sources. 
. Results. It is shown that the expected interstellar magnetic fields substantially overpredict even these 7-ray upper limits. Only rather strongly 
' amplified magnetic fields could be compatible with such low 7-ray fluxes. However this would require a strong component of accelerated 
nuclear particles whose energy density substantially exceeds that of the synchrotron electrons, compatible with existing theoretical acceleration 
^\ ' models for nuclear particles and electrons. 

' Conclusions. Even though the quantitative arguments are simplistic, they appear to eliminate simplistic phenomenological claims in favor of 
00 , a inverse Compton 7-ray scenario for these sources. 

q 
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1. Introduction the field fluctuations are excited by the accelerating particles 

^, . , , , ,. , themselves and since the pressure of these particles (typically 

The question, whether the very high energy (VHE) (h^ > , , ... , n ■ u 1 4.1, 

^ . \ / V 7 comparable with the thermal pressure) is reacting back on the 

100 GeV) 7-ray emission of the Galactic supernova remnants ^, , , • . 1 1 j . i. 

' ^ thermal plasma, this strongly coupled system becomes a com- 

(SNRs) implies a sufficiently large production of nuclear cos- , ^ ^■ j • . 1 r ^i. i, j 

J o ±- p^Qx problem of nonlinear dynamics, not only tor the charged 

mic rays (CRs) - of the same order as that required to replenish . , . 1. . 1 r 1 . . .c u j 

-1 ±- particle components but also lor the electromagnetic field and 

the Galactic CRs - IS one of the key problems addressed by 7- r, . j 1 . . • 1.1 r 

•' ^ ' Its fluctuations. Models suggest that a sizeable fraction of the 

ray astronomy. There are two ways to deal with this question u j j • 1 • ^ c j • » 

■' ■' ^ entire hydrodynamic explosion energy will be transformed into 

in the investigation of an individual SNR. . . , t,.- . . .u cxrr. 

energetic particle energy. This suggests that the SNRs are in- 

The first approach is a theoretical one. It uses a nonUnear ^^^^ ^^^^^^^ ^j^^ ^^j^^^.^ 

combination of gas dynamics (or eventually magnetohydrody- 

namics) for the thermal gas/plasma and kinetic transport theory nuclear charged particles and electrons can be acceler- 
for the colUsionless, nonthermal relativistic particle component ^^^d to achieve nonthermal momentum distributions. The ener- 
that is coupled with the plasma physics of the electromagnetic S^tic electrons show their presence through synchrotron emis- 
field fluctuations which scatter these particles. In the environ- ^ion from radio frequencies to hard X-ray energies. They may 
ment of the coUisionless shock wave of a supernova explosion ^l^o interact with diffuse interstellar radiation field photons, 
this allows the description of diffusive shock acceleration of the l^^e the Cosmic Microwave Background, to produce high en- 
energetic particles which are originally extracted from the ther- ^rgy 7-rays in inverse Compton (IC) collisions. The injection 
mal gas and thus injected into the acceleration process. Since of electrons into the acceleration process is however poorly un- 
derstood quantitatively. Even assuming that the electron mo- 
Send offprint requests to: H.J. Volk mentum distribution at high particle energies is only propor- 
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tional to the total number of electrons injected per unit time 
at the shock, the amplitude factor of the electron momentum 
distribution is therefore not known from theory. It is typically 
inferred from the measured synchrotron spectrum produced by 
the accelerated electrons by assuming a mean strength of the 
magnetic field. This will be a key question in the discussion of 
this paper. 

The injection of nuclear particles from the suprathermal tail 
of the momentum distribution produced in the dissipative shock 
transition is much better understood, because the same mech- 
anism that produces scattering fluctuations for higher-energy 
nuclei - essentially a beam instability from the accelerated nu- 
clei that is so strong that asymptotically the particles scatter 
along the mean field direction already after one gyro-period - 
also works at injectio n energies. In addi tion, where heavy ion 



injection takes place (IVolk et al.L 120031) . it is to be expected 



that these ions dominate the nonthermal energy density behind 
a strong shock like in a young SNR. Then the main nonlin- 
ear shock modification consists in a weakening of the quasi- 
discontinuous part of the shock structure, associated with a 
broad shock precursor Low energy particles - ions and elec- 
trons - in the accelerated spectrum are then only accelerated 
at this weaker subshock and this implies a significantly softer 
momentum spectrum at low energies than at high energies. This 
physical effect is visible in the radio part of the electron syn- 
chrotron spectrum and therefore a quantitative indication of the 
degree of shock modification. It provides a means to determine 
the injection rate of nuclear particles, de facto of protons, from 
the radio synchrotron spectrum. In all cases, where the syn- 
chrotron spectrum of SNRs was measured, this softening was 
observed. Together with the nonlinear theory of acceleration, 
and in the strong scattering limit, this determines the nonther- 
mal pressure P^, which turns out to be comparable with the ki- 
netic pressure pV^ of the gas. Here Vg and p denote the shock 
velocity and the the upstream mass density, respectively. Using 
therefore the synchrotron measurement, the nonthermal quan- 
tities can be determined from theory. The exception is at first 
sight the mean magnetic field strength. However, it needs to be 
consistent with the overall form of the synchrotron spectrum, 
from radio to X-rays, and with the X-ray synchrotron morphol- 
ogy that depends on the effective strength of the magnetic field. 
In this way an interior effective magnetic field strength is de- 
termined. It is typically an order of magnitude larger than the 
MHD-compressed upstream field strength. This amplification 
of the magnetic field is a characteristic of the effective acceler- 
ation of CR nuclei in a SNR, because it can only be the result 
of strong acceleration of nuclear particles. The pressure of the 
accelerated electrons - also for the cases discussed below - is 
more than two orders of magnitude below pV^ . 

The question is then, whether the observed 7-ray emission 
is also dominated by nuclear particles through their inelastic, 
7r° - producing collisions with thermal gas nuclei. This need 
not be the case if the target density of the thermal gas is very 
low, despite the fact that the energy density of the accelerated 
nuclear particle component is very high, in fact comparable to 
the thermal energy density. 

Using this theoretical approach (for reviews, s e e e.g 



the investigation of half a dozen of young Galactic SNRs has 
shown that the nuclear CR production is in all cases so high that 
the Galactic SNRs are viable candidates for the Galactic CR 
population up to particle ener gies ~ 10^^ eV, well abo ve the 



so-called knee in the spectrum ( IBerezhko & VoUd 120070 . Even 



though important details are open to debate because the time- 
dependent ev olution of a point explo sion can only be calculated 
numerically ( Berezhko et al. . 19961) . we believe that this result 



is quite a robust one. 

However, from a strictly observational point of view, the 
hadronic nature of most of the SNR 7-ray sources is not proven 
this way. This might ultimately be possible in a direct way 
with a very sensitive neutrino detector. The remaining ques- 
tion whether the Galactic CR population has a SNR origin then 
still requires the consistency of the observational result and the 
theoretical picture. 

As far as 7-ray observations are concerned, there is also 
a different approach, basically phenomenological. It consid- 
ers the question, whether and to which extent the hadronic 
or leptonic origin of the measured 7-ray emission can be 
decided by favoring either one mechanism at the expense 
of the other directly from the data. For example, it can ask 
the question whether the necessarily limited dynamical range 
of the observed 7-ray emission allows a distinction between 
a hadronic and a leptonic scenario. Or it can ask whether 
observations in other wavelength ranges tend to empirically 
contradict the theoretically favored scenario of a predom- 
inantly nuclear energetic particle energy density. A possi- 
ble topic consists in the interpretation of spatial correlations 
in res olved 7-ray SNRs, lik e those noted in RX J1713.7- 
3946 jAharonian et all l2006l) and RX J0852.0-4622 (Vela Jr.) 



jAharonian et al. 

L l2007h . The correlation of the hard X-ray syn- 



chrotron emission with the VHE 7-ray emission features might 
be considered to favor energetic electrons to produce both 
emissions. Discussions of the above and similar issues have 



recentlv been given for instance in Aharonian et alj (20061); 


Porter et al. (l2007l); lAharonian et alj (2007b; 


Katz & Waxman 


(2008); Plagal (2008b. and Berezhko & VoUd ( 


2008b. However, 



the complexity of the configurations that characterize these ex- 
tended sources introduces severe uncertainties. They arise from 
the poorly known structure of the circumstellar medium, which 
could be partly due to the strong winds expected from the 
progenitor stars or could be partly pre-existing in the form of 
neighboring interstellar clouds, affected by the progenitor and 
its subsequent explosion. 

We shall add in this paper such a phenomenological argu- 
ment. It concerns the spatially integrated synchrotron emission 
spectrum for the simplest available objects, the remnants of 
the three young type la SNe, observed in VHE 7-rays. Even 
though only upper limits exist from the H EGRA, H.E.S.S . 
and C ANGAROO exp eriments for SN 1 ( )06 d Aharonian et al 



2005b, Tycho's SNR dAharonian et all I2OO1]). and Kepler's 
SNR dEnomoto et al.L I2OO8I: lAharonian et all l2008b . they can 



Malkov & Drurvl I2OOII; IvHm 120041; Berezhkol I200I l2008i) 



nevertheless be used to estimate lower limits to the effective 
mean magnetic field strengths in the SNR that are consis- 
tent with the observed spatially-integrated synchrotron spec- 
tra. These somewhat naively estimated magnetic fields are then 
compared to the expectations for these types of SN explosions. 
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The large discrepancies found disfavor leptonic scenarios for 
these objects. 

2. Simple synchrotron and IC modeling of the 
integrated emission 

The expected synchrotron spectral energy density (SED) at 
distance d frorn a SN R is given by the expression (e.g. 
Berezinskii et al. 1 19901) 



,f^^syn 3x10- 



dE 



d^rB± 



X / dpp foir,p)g . 



E 



(1) 



in erg/(cm^s), where 



g{y) K^/^{y')dy' , 

Kf^{y) is the modified Bessel function, E is the photon energy, 
= 3eB±p'^ /[ATT{mcc)^], and B± is the interior magnetic 
field component perpendicular to the line of sight. 

We shall use here an approximation that averages over 
the line of sight directions. A precise analytical integra- 
tion involving Whittacker's function has been given by 
Crusius & Schlickeiser ( 19861) which is in turn closely approx- 
imated by substituting 

(2) 

into Eq.([T]). Here is the strength of the interior field which 
results from the MHD-compression of the upstream magnetic 
field Bo and subsequent de-compression in the interior (see 
below). The strength of Bq is denoted as Bq. 

The spatial integral in Eq.([T]) extends over the volume V 
of the SNR, as given by the observed synchrotron morphology, 
and the calculated synchrotron SED has to be compared with 
the observed SED. 

Our starting point for a simplified model is the assumption 
that B± in the form of Eq.© can be taken as a weighted mean 
value \/2/3{Bd) out of the spatial integral of Eq.([T]i. The post- 
shock value of B^/Bq is locally between 1 and a, where cr > 1 
denotes the overall shock compression ratio. Since the interior 
field strength is lower than the postshock field strength, we have 
for this mean interior field strength: (Bd) < o-Bq. 

If we investigate the possibility that the accelerated parti- 
cles are electrons alone - implying a purely leptonic origin of 
the VHE 7-ray emission - then we have to consider a test parti- 
cle problem with cr = 4. In the same sense Bq should be equal 
to the strength of the interstellar magnetic field, i.e. equal to 
a few yuG. Values of Bq = 3/iG and Bq = 5/iG then imply 
3 < (i?d) < 12/iG and 5 < (i?d) < 20juG, respectively. 

As a second approximation we shall also assume that the 
volume integral of fc{r,p) equals the product of V and an elec- 
tron distribution 



d^rfc{r,p) = VAp " x cxp(-p/pniax), 



in the form of a power law with an exponential cutoff at Pmax- 
The index a = 4 again corresponds to a test particle spectrunu. 

In this sense the two parameters A and Pmax can be ap- 
proximately fitted from the known radio and X-ray synchrotron 
data as a function of Bq. In fact, because of the exponential be- 
haviour of the cut-off of the electron momentum distribution 
the only sensitive parameter turns out to be A. The fact that the 
observed radio synchrotron spectra are softer than implied by 
a distribution with a = 4 suggests that the pure electron accel- 
eration model, for the sake of argument considered here, is not 
the physically correct model. However, a pure electron model 
is necessarily one with a = 4, even if it does not optimally fit 
the form of the observed synchrotron spectrum, but rather only 
its amplitude. 

Next we calculate the IC SED from these same electrons 
in the 2.7 K cosmic microwave background (CMB). This 7-ray 
SED can be written in the form: 



E^ — 2_ 

dE 



I d^r I deuphie) 



X / dpp^a{ee,E,e)fe{r,p) 

"'Pmin 

in erg/(cm^s), where dBlumenthal &Gouldl[l970l) 

3fTT(TOeC^)^ 



(4) 



4ee2 



2qlliq+{l + 2q){l-q) + 0.5 



{rq)^l-q) 



1 + Tq 



(5) 



is the differential cross section for the up-scattering of a photon 
with incident energy e to energy E by the elastic collision with 
an electron of energy ec. 



1 



(6) 



is the blackbody spectrum of the CMB, h ~ 2TTh and ks are 
the Planck and Boltzmann constants, respectively, T = 2.7 K, 
(It = 6.65 X 10^2^ cm^ is the Thomson cross-section, q = 
E/[r{ec — E)], r = 4eee/(moc2)^ , and pmin is the minimal 
momentum of the electrons, whose energy ec is determined by 
the condition q = 1. 

We neglect here nonthermal Bremsstrahlung emission 
which turns out to be unimportant for all the cases considered 
below. 

Since the CMB is uniform, we can without further approx- 
imation use Eq.(l3]i to express the 7-ray SED in terms of the 
parameters A and Pmax- The results are given in Fig.[T^ for SN 
1006, in Fig. [TJ) for Tycho's and in Fig. [T]; for Kepler's SNRs 
for various values of (B^). 



(3) 



' For internal magnetic field strengths in excess of 100/iG such 
a model distribution would have to include a high-energy part 
of the spectrum that is softened by synchrotron losses, see e.g. 
iBerezhko et"ai] j2002l) . However in the present context, that does by 
assumption not include massive nuclear particle acceleration, such 
field strengths are not expected to occur. 
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Fig. 1. The overall (spatially integrated) nonthermal spectral energy distribution (SED) as a function of photon energy E. The 
lower-energy part shows the simple fit to the observed synchrotron-SED, cf. Eq.(|3]i, for various values of the internal field strength 
B± in fj,G, cf. Eq.(|2]i. The synchrotron fit is essentially independent of mean strength (Bd) of the internal field. The high-energy 
curves show the inverse Compton-SED in the CMB for the various field st rengths, (a) for SN 1006: The blu e radio data are from 
Revnolds d 19961) whereas the ^ reen and red X-ray data are from RXTE jAUen et all Il999 l). and Suzaku dBamba et al.l I2OO8I) . 
respectively. The Chandra data ( lAllen et all 120041) are very similar to the Suzaku data and can be treated as indistinguishable in 
the present context. Also given are the upper limits from H.E.S.S.dAharonian et al.ui2005i) and EGRET (^Naito et al.uil999i) . (b) 
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3. Discussion 

It is clear from the outset that the approximate nature of the 
models used limits the impact of the conclusions to be drawn 
from these results. On the other hand, most of the arguments 
that have been used in the past regarding the alternative be- 
tween a hadronic and a leptonic interpretation of VHE 7-ray 
results have used such one box approximations. The only alter- 
native would be full time-dependent solutions of the governing 
system of equations, discussed in the Introduction. 

However, with this proviso, the results are surprisingly 
clear. For all three sources magnetic field strengths Bo lower 
or equal to the expected interstellar magnetic fields of 3 — 5yuG 
substantially overpredict even the existing 7-ray upper limits. 

For a shock that excites MHD fluctuations only weakly if 
at all, because of the assumed lack of acceleration of nuclear 
particles, the interior gas flow will be essentially laminar and 
adiabatic. Taking into account that in such an approximately 
laminar gas flow the minimum strength of the internal mag- 
netic field will always be lower than the strength of the up- 
stream field and that the weighted average field strength av- 
erage field strength is considerably lower than the maximum 
field strength (over the quasi-circular shock surface) immedi- 
ately behind the shock, a more realistic estimate for the value 
of (Sd) would be to put a ^ 1. This would imply that the 
(i?d) - values, given in FiglU roughly equal the values of Bq. 
This means that already the curves for B^ ~ lOyuG in the fig- 
ures assume an unrealistically high ambient interstellar field 
strength, larger than the interstellar average. Yet they overpre- 
dict the IC 7-ray flux by at least one order of magnitude in 
comparison with the observed total 7-ra y upper limi t alrea dy 



for the very low-density object SN 1006 (lAcero et all l2007n 



that would therefore be expected to be located in a lower than 
average interstellar magnetic field as well - and by much more 
for the two other sources. 

Existing theoretical solutions for the overall particle accel- 
eration in these three sources take into account the amplifi- 
cation of the magnetic field by the accelerating nuclear par- 
ticles whose energy density becomes comparable to the ki- 
netic energy of t he incoming gas flow, as seen in the frame 



of the shock (e . g. iKsenofontov et al.L 120051: IVolk et al.L 12005 



Berezhko et alil2006l) . Only then it seems possible to not over- 



predict the leptonic flux. At the same time the 7-ray flux is 
dominated by the hadronic flux, even though in SN 1006 only 
by a small margin. 

We note here that besides nonhnear amplification due to 
CRs the magnetic field in SNRs can also by amplified by other 
mechanisms. These are Rayleigh-Tay lor instabilities at the con- 
tact discontinu ity between the ejecta and the shocked circum- 
stellar gas (e.g. lWang & Chevalien . l200lh and the vorticity gen- 
eration that results from the shock running into possibly pre- 
existing density inhomogen eities of the circumstellar medium 
(iGiacalone & Jokipiill2007l) . However, the common feature of 
these mechanisms is that they act only in the downstream re- 
gion and produce their main effect at a substantial distance 
behind the shock, that is outside the CR acceleration region. 
Therefore these mechanisms do not influence the CR acceler- 
ation process, in particular the maximal particle energy. Since 



accelerated CRs in young SNRs are concentrated in a thin layer 
near the shock front, these two mechanisms also hardly influ- 
ence the properties of the nonthermal emission, produced by 
CRs, except possibly that of the highest-energy CRs. 

4. Conclusions 

Simple one box approximations indicate that a leptonic sce- 
nario for the 7-ray emission from the three known Galactic 
type la SNRs SN 1006, Tycho's SNR and Kepler's SNR sig- 
nificantly overpredicts the 7-ray flux, even when compared 
to the existing upper limits from observations. The calculation 
makes direct use of the observed synchrotron emission spectra. 
Even though the arguments are simplistic, they appear to elim- 
inate equally simplistic phenomenological arguments in favor 
of such a scenario. Any positive argument in favor of a purely 
leptonic scenario would therefore have to be based on a full so- 
lution of the governing nonlinear equations. From our results, 
however, we believe that such a positive argument can not be 
made. 
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